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ABSTRACT 



The infrared vibrational spectra of HOD and DOD dissolved in carbon 
tetrachloride and carbon disulfide were obtained and compared with the 
same spectra of HOH. The spectra of these three water species were com- 
pared to those of the vapors and the similarities noted. The bending 
mode, 2 > °f DOD in carbon disulfide was observed. Temperature effects 
on the spectra were also studied. To analyze the data a simple theoret- 
ical vibrational treatment was used providing a computation readily 
adaptable to digital computers. 
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1. Introduction 



Infrared Spectra of HOH 

The infrared spectrum of water has been intensively studied in 
vapor phase, liquid phase, and in dilute solutions of non-polar solvents. 
The vapor phase spectrum has been analyzed by many investigators and in 
considerable detail first by Darling and Dennison, j^lj 

The water molecule is well known to have the form of an isosceles 
triangle with equivalent OH bonds, an example of symmetry. Such a 

molecule possesses three normal modes of vibration, all of which should 
be active in both the infrared and Raman spectra. The H-O-H bond angle 
is reported to be 104.5 ± 0.10° with an equilibrium bond distance of 
0.9580A. jl-3j The bond angle of DOD is known to be identical. j\j 
Benedict, Gailar and Plyler present the most recent definitive 
infrared spectral data for HOH, giving the values of the fundamental 
water vapor vibrational frequencies as: y ^ = 3656.65 cm -1 , y 2 ~ 1594.56 
cm“l, and = 3755.79 cm“\ j\J These vibrational modes may be de- 
picted with reasonable accuracy as shown in Figure A. 
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Figure A 

Herzberg summarizes previous values of infrared and Raman data for 
the various water molecules HOH, HOD, and DOD. 

Liquid water exhibits the characteristic very broad, intense infra- 
red absorption found for hydroxy containing compounds such as alcohols 
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and carboxylic acids. This band is attributed to associated or hydrogen 
bonded species. 1 ^ 6 j 

Infrared Spectra of Water in Non-Polar Solvents 

Upon dilution of the OH containing compounds with non-polar solvents 
such as carbon tetrachloride or carbon disulfide, a new, narrow absorp- 
tion band appears at a higher frequency than that of the original associa- 
ted species. |7 -9j This new absorption is assigned to the stretching 
frequency of the unassociated OH bond. 

Infrared spectra of dilute solutions of HOH in carbon tetrachloride 
were studied by Borst et al., who observed two absorptions in the OH 
stretching region at 3702 cm"'* - and 3611 cm - -*-, jioj These investigators 
suggested that the 3702 cm"'*- absorption be assigned to 9^3 and the 3611 
cm - -*- band to *V ^ . Later Fox and Martin confirmed these observations at 
frequencies of 3705 cm"-*- and 3614 cnT* and stated that these absorptions 
were, respectively, the 9^3 and 9 ^ vibrations of the monomeric HOH 
species in carbon tetrachloride solution. Jllj Reisinger has confirmed 1 
these observations assigning 3705 cm"^ to ^3 and 3613 cm - -*- to 0^. jl 2j 
None of these investigators reported an observation of the *y ^ or bend- 
ing frequency. This is to be expected, since in the path lengths used 
for these investigations carbon tetrachloride and carbon disulfide 
absorb intensely in the region near 1595 cm - ^, the vapor phase 9^2 fre- 
quency. 

To avoid these solvent absorptions the present investigation utilized 
DOD and HOD whose gas phase V 2 frequencies are 1178 cm - -*- and 1402 cm - -* 
respectively. This yielded an interference free spectrum only for the 
I /2 mode of DOD in carbon disulfide, the 9 / 2 of HOD being obscured in 



all solvents used. 



Fox and Martin were able to separate graphically a third band in 
the HOH spectrum centered near 3810 cm"*, llj Noting that this absorp- 
tion appeared to follow the general shape of rotation lines in this 
region they proposed that this absorption was the envelope of the rota- 
tion lines, Reisinger also noted this third band placing it at 3804 
cm" 1 . 12 Additionally, he separated a fourth absorption in this region 
centered at 3540 cm" 1 . These third and fourth absorptions were also 
noted by Reisinger in an investigation of a dilute solution of HOH in 
carbon disulfide. |12j The nature of the absorptions appeared identical 
to those of HOH in carbon tetrachloride except for a slight shift of 
the frequencies to 3796 cm" 1 and 3547 cm" 1 . Earlier investigators do 
not indicate that their studies covered the low frequency region (near 
3500 cm" 1 ) so that this fourth band (3540 cm" 1 ) was unobserved by them. Jllj 
The 0^2 frequency for HOH in carbon tetrachloride and carbon disul- 
fide would be expected at about 1595 cm" 1 (the vapor phase since, 

as is reported later, the 9^2 f re( 3 uenc y for DOD in carbon disulfide is 
very little shifted from its vapor phase value of 1178 cm'^. However, 
ry 2 has not been observed for HOH in carbon tetrachloride or carbon di- 
sulfide. 

Temperature Effects 

Early investigators of the spectra of OH bearing species noted 
that the relative intensities of the OH stretching absorptions were 
functions of temperature. They attributed the effect to temperature 
induced changes in monomer-polymer equilibrium or to solvent effects, 
i.e., a change in the molar absorption coefficient. ^8,13 

More recently, however, Slowinski and Claver have noted that for 
pure liquids the intensities of OH stretching absorptions as well as 
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C-H, C-Cl and other bond stretching absorptions depend markedly on temper- 
ature. jl4j These workers show clearly that both the maximum optical 
densities and the integrated intensities of these stretching absorptions 
decrease with rising temperature. Their work was rapidly followed by 
the efforts of other investigators who have extended and refined these 
observations, studying hydrogen bonded equilibria in dilute solutions. 
jl5,16j Hughes, Martin, and Coggeshall demonstrate that the intensity 
of the absorbance of the unassociated OH stretching bond decreases lin- 
early with increasing temperature, jisj 

Calculations 

The theoretical bases for correlating vibrational frequencies of 
molecules with their parameters (masses and geometries) has been known 
for many years. Wilson, Decius, and Cross have thoroughly delineated 
the methods for these calculations. jl7j The matrix methods presented 
by these authors are readily adapted to digital computer technology 
and presented an opportunity to determine if there was a simple theoret- 
ical vibrational treatment consistent with the results of this investi- 
gation. In one simple treatment of a non-linear triatomic molecule 
such as water, the three fundamental vibrational frequencies are deter- 
mined by four force constants and the molecular parameters. These three 
frequencies are, of course, insufficient information for solution of 
the equations determining the four force constants. However, applica- 
tion of the Born-Oppenheimer approximation permits the assumption that 
the potential function is unchanged by isotopic substitution. This 
allows the six fundamental vibrational frequencies of HOD and DOD to be 
used in determining these same force constants. In the two calculations 
made here the first used the five observed frequencies of HOH and DOD. 



The second 
a total of 



used, in addition, the two observed frequencies of HOD for 
seven. 
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2. Experimental 

Infrared Absorption of the Water Molecules HOH, DOD, and HOD 

The Perkin-Elmer Model 221 recording spectrometer with a sodium 
chloride prism-grating combination was used in the double beam mode of 
operation. The spectrometer was calibrated in the regions of interest 
using the vapor phase spectra of HC1, HBr, NH3, C0 2 , and HOH as indicated 
in Table I. The pressures were chosen for the purpose of direct compar- 
ison with the reference values given by the Commission on Molecular 
Spectroscopy, jl 8^ All gases used except the HOH were Matheson Co. bot- 
tled gases. The HOH spectrum was obtained from atmospheric vapor using 
the spectrometer in single beam operation. The standard deviations from 
the mean of the observed frequencies are given in Table I for each range. 
A mean correction was determined for each range and applied to the ob- 
servations to yield the reported values. 

Connecticut Instrument Corporation 50 mm path length glass cells 
were used with BaF^ windows. These cells have a capacity of approxi- 
mately 30 ml. Only glass, Teflon, and BaF 2 were in contact with the 
sample in these cells. 

The carbon tetrachloride and carbon disulfide used as solvents 
were Fisher Scientific Co. spectral grade reagents. No additional puri- 
fication was necessary. It was necessary to dry the solvents as com- 
pletely as possible using the disappearance of the 9^3 and 7^ HOH fre- 
quencies as a minimum criterion of dryness. Previous investigators 
have used a variety of techniques for drying, most of which proved un- 
satisfactory when used in this laboratory. ^9,12,19 These included 
drying over CaCl 2 , Na 2 S0 3 , Na 2 S0^, bubbling with dry air, and distilla- 
tion. P2°5 was a successful drying agent, but its use involves cons id- 
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TABLE I 



Calibration Data 



Standard 


Frequency (cm"* 1 ) 


(cm" 1 ) 


Pressure (mm) 


HCl 


2680-2900 


0.5 


120 


HBr 


2600-2720 


0.5 


140 


nh 3 


1160-1200 


0.3 


50 


h 2 o 


3745 


0.9 


* 


co 2 


3620 


0.9 


400 



* Atmospheric vapor 

erable tedious handling. Linde Co. type 4A molecular sieve in bead 
form was used in this investigation, the beads being placed directly 
in a glass stoppered bottle of the solvent. In order to prevent leak- 
age of moisture into the sample, the bottle was placed in a dessicator 
containing P 2 O 5 . No clian g e was observed in the dried solvents or the 
molecular sieve. The solvent spectra were also unchanged over the re- 
gion 4000 to 750 cm -1 except, of course, for the HOH bands. All trans- 
fers were made using a syringe to minimize contact with atmospheric 
moisture . 

Dry carbon tetrachloride with a small amount of added DOD produced 
an observable HOD absorption in as short a time as fifteen minutes when 
placed in the glass-stoppered BaF 2 cell. However, dry carbon tetra- 
chloride alone did not produce an HOH absorption for a minimum of several 
times as long. It was proposed that exchange of D for H with the glass 
cells may have been responsible for this phenomenon. To prevent con- 
tamination by absorption of atmospheric moisture a few granules of 
molecular sieve were placed in the dry reference cell sample. This 
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effectively prevented contamination during a span of several hours, which 
was the longest period observed. All spectra were recorded with an 
identical cell filled with dry solvent inserted in the reference beam 
of the spectrometer. 

Saturated solutions of HOH and DOD in both solvents were prepared 
by adding excess HOH or DOD to the dried solvent and allowing a minimum 
of 24 hours for equilibration. The DOD was obtained from Bio-Rad Lab- 
oratories and stated to be 99.83 mole percent DOD. 

The various sample solutions were prepared as a given percent of 
saturation by dilution of the saturated solution with dry solvent. 

The equilibrium for the reaction 



HOH + DOD * 2 HOD with 



[HOD] 2 

[hohj qx)d] 



is known to be very rapid. The equilibrium constant has been determined 
as K = 3.96 at 25°C. J2oj For a given sample we have 

[hod] + [hoh] + [dod] + C 

where C is a constant. Noting that 



[hod] 2 = k [hoh] [dod] 

it can be seen that the condition for HOD to be a maximum is 

[hoh) = [dod] . 



So then 



[hod) 



max 



rr 

2 + VK' 1 " 



and for K = 3.96 the maximum [HOD] is about 50 mole percent, which is 
a ratio of 1:2:1 for the various water species. This maximum HOD con- 
centration exists for an equimolal mixture of HOH and DOD. The values 
in Table II give the observed frequencies for the HOD spectrum in car- 
bon tetrachloride. The notation of Benedict et al. has been adopted 



TABLE II 



Frequencies of vibrational absorptions of various species 
of water (in cm" 1 ) 



Vapor Phase 





HOH 


DOD 


HOD 


*1 


3657 


2671 


2727 


^2 


1595 


1178 


1402 


^3 


3756 


2788 


3707 






Carbon Tetrachloride Solutions 






HOH 


DOD 


HOD 


« — \ 
V 


3614 


2642 


2692 


v 2 


- 


- 


- 


v 3 


3705 


2751 


3660 






Carbon Disulfide Solution 






HOH 


DOD 






3603 


2635 




^2 


- 


1180 




?3 


3695 


2744 





Values given for vapor phase spectra are those of Bene- 
dict et* al. (rounded to the nearest cm" 1 ). 4| 
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in describing the HOD vibrational bands, namely = OD stretching, 

0^2 = bending, and ^3 = OH stretching vibrations, 

There is considerable overlap of the y^ HOD frequency, 2692 cm” 1 , 
with those of the DOD spectrum and of the HOD frequency, 3660 cm” 1 , 
with the HOH spectrum. Consequently, the optimum concentration of HOD 
is one that will minimize the interfering bands while allowing sufficient 
intensity of the HOD absorption. A 9:1 misture of HOH and DOD will 
yield approximately an 81:18:1 ratio of H0H:H0D:D0D. This 18:1 ratio 
of HOD to an interfering species is an order of magnitude improvement 
over the previous ratio of 2:1 at maximum mole percent HOD. 

Temperature Effects 

In order to observe the effect of temperature on the spectra 
special cells were constructed. These were 50 mm path length glass 
cells with a coolant jacket as an integral part of the cell. Water 
from a constant temperature bath was used as the regulating fluid. 

Silver chloride windows were held in place on the ends of the cell body 
by aluminum discs fastened to each other with screws. Internal cell 
temperatures were measured with a thermocouple. Carbon tetrachloride 
solutions were studied over a range from 3°C to 60°C. Due to the high 
vapor pressure encountered in this temperature range, carbon disulfide 



solutions were not studied. 



FREQUENCY (CM ) 
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3. Results and Discussion 



The results of earlier workers were confirmed upon analysis of the 
spectrum of a dilute solution of HOH in carbon tetrachloride. |l0-12j 
This is shown in Figure 1. Two of the fundamental vibrational frequen- 
cies were identified as 9^ at 3614 cm” 1 and 9^3 at 3705 cm” 1 in close 
agreement with previously reported values. ^1 1 , 1 2 J Easily identified 

also were the broad sidebands centered at 3804 cm -1 and 3540 cm - " 1 . These 
sidebands were graphically separated from the spectrum. They were also 
observed by this technique by Reisinger in a solution of HOH in carbon 
disulfide with centers at 3796 cm” 1 and 3547 cm"" 1 respectively. jl2j Fox 
and Martin observed only the 3804 cm -1 sideband considering it to be the 
envelope of rotation lines of the water molecule. JllJ Several possible 
sources for the 3540 cm” 1 band were considered by Reisinger including: 

(1) the possibility of Fermi resonance between 2 and (2) that it 

is 2*V 2> an ^ (3) that it is the envelope of rotation lines. Jl2j 

Fermi resonance, a quantum mechanical phenomenon, is the result of 
interaction of closely spaced energy levels. In this process these 
levels are spread apart and associated transitions tend to become more 
nearly equal in intensity. A necessary condition for Fermi resonance 
is that the levels involved be of the same symmetry species or share 
at least one species in common. Since 2 0^2 and 0 /\ share the fully 
symmetric species , Fermi resonance is possible if the energy levels 
are sufficiently near one another. 

A determination of the absorption of D0D in carbon tetrachloride 
at varying percentages of saturation at 28oc confirmed that the bands 
at 2751 cm” 1 and 2642 cm” 1 should be assigned to the and fre- 

quencies of the D0D molecule. These results are shown in Figure 2. It 



is observed that the structure of the DOD absorption is completely anal- 
ogous to that for HOH including the two low intensity sidebands. The 
spectrum of HOD in carbon tetrachloride shown in Figure 3 for an eighteen 
to one HOD to DOD mixture is also observed to have the same basic 
structure for Y^ and Yy 

Since the solvents investigated are opaque in the region of the 
HOH bending vibration near 1595 cm'^, a search was made for the bending 
frequencies of DOD and HOD. The HOD bending mode was found to be ob- 
scured in 50 mm path lengths of both carbon tetrachloride and carbon 
disulfide. For DOD in carbon disulfide it was possible to observe a 
low frequency band at 1180 cm“l which was assigned as the bending vibra- 
tion, 7^2* In view of this appearance of 9^2 ver Y close to the gas 
phase frequency of 1178 cm'^, it appears that the bending vibration of 
the water molecule is almost entirely unaffected by solution in carbon 
tetrachloride and probably therefore in other non-polar solvents. 

Hornig, using infrared data obtained from systems of dilute HOD 
in both liquid HOH and DOD, states that the shape of the HOH absorp- 
tions in the stretching regions results from overlap of Y\> an< ^ 

2 V 2 , and that the intensity variations are caused by Fermi resonance 
between 2 Y 2 and [21,22j This type of system exhibits a decrease 

of about 300 cm - -*- in , and an increase of about 50 cnT^ in Y 2 w ^ en 
compared to gas phase frequencies. These shifts bring 2 7^, and 7 ^ 
sufficiently near one another to allow the possibility of Fermi resonance. 

Considering the present case of HOH in carbon tetrachloride, it 
seems by analogy with DOD in carbon disulfide, that the shift of ^ 
from its gas phase frequency is essentially zero. The shift of is 
a decrease of about 50 cm"^. The 2 mode for HOH vapor is observed 
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at 3151 cm \ lowered from the harmonic potential value of 3190 cm”^. 
Considering the unshifted position of V 2 non-polar solvents it is 

reasonable to conclude that it is very unlikely that 2 V 2 will exceed 
3190 cm"'*' in solution in non-polar solvents. This would eliminate 2 7^ 
as an explanation of the 3540 cm“^ sideband. These considerations lead 
to a minimum spacing between 2 2 and of about 400 cm*" 1 implying 

that Fermi resonance is unlikely as an explanation of the broad, low 
intensity band at 3540 cm”-*-. This disposition of 2 7^2 and Fermi 
resonance leaves as the remaining plausible explanation of the band 
the possibility that it is an envelope of rotational lines. 

The persist ance of rotational structure of a vibrational band of 
a molecule passing from a vapor to a liquid phase has been given little 
theoretical treatment. Longuet -Higgins expects that hindered rotation 
may persist through this transition, and that for small, nearly spheri- 
cal molecules the rotational envelope should be recognizable in favor- 
able cases. j^2 3 j The interpretation of the 3804 cm~l and 3540 cirT-*- side- 
bands as rotational envelopes is in this view a reasonable assignment. 

If this assignment is correct it might be possible to observe a 
change in these envelopes as a function of temperature. The spectrum 
of an asymmetric rotor such as water will not show the simple changes 
observed for a diatomic molecule. However, an increase in temperature 
giving an increase in the population of higher energy levels would pro- 
duce the same basic change in a gross feature such as the rotation en- 
velope; an extension of the envelope to a wider frequency range with 
greater relative intensity for the extremes of the envelope. In con- 
trast to this expected increase of intensity with temperature, the re- 
sults of Hughes et al. (that the intensity of the OH absorbance decreases 



with increasing temperature) leads to the possibility of either an in- 



crease or decrease. In fact, no change in intensity of the envelope 



was was identified over the 3°C to 60°C range observed. Some tentative 
changes were noted but were not reproducible. The intensities of the 
vibrational absorption maxima did decrease with increasing temperature 
in approximately the manner described by Hughes for all three species 



of water molecules in both solvents. Due to the difficulty of 



accurately measuring integrated intensities, only maximum optical den- 
sities were noted. Using the cells described, maintenance of a uniform 
temperature was not possible. It would be desirable to use an apparatus 
of the type used by Hughes: a specially constructed metal compartment 

which maintains an accurate, uniform temperature over a wide range. Jl5J 
Some manner of sealing the system under pressure would also be useful. 

The fact that the 3804 cm** 1 and 3540 cm~l sidebands did not show 
a detectable decrease of intensity over this temperature range indicates 
that these bands are very probably not overtone or combination bands 
since the temperature effect would have markedly reduced their intensi- 
t ies . 

It seems, then, that although there is little direct evidence that 
these two sidebands are rotational envelopes, the dismissal of all other 
proposed sources for the bands leaves this as the most probably inter- 
pretation. 

The hygroscopic natures of carbon tetrachloride and carbon disul- 
fide are exhibited in the great difficulty encountered in drying the 
solvents and in the rapid re-absorption of the water. A glass stoppered 
cell of carbon tetrachloride began to show HOH absorptions in periods 
on the order of one hour. Almost paradoxically, however, the solvent is 
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The tenacious grip 



saturated by a mere one part in 10 ^ of water. |19 
of carbon tetrachloride and carbon disulfide on water may well be ex- 
plained by interstitial trapping. Such a mechanism would reconcile the 
observed difficulty in removing trace amounts of water with the low 
solubility while also allowing the water molecule to remain essentially 
free of interactions. 

A normal coordinate analysis of the data was performed using the 
methods of Wilson, Decius, and Cross, jl 7^ The experimental values of 
the fundamental vibration frequencies were used to calculate the quad- 
ratic force constants in the potential energy expression. The constants 
treated are those for a modified valence force treatment. In the case 
of water, there are four constants. These are for OH bond stretching, 
H-O-H bond angle deformation, and two constants for bond interactions 
(stretch-stretch and bend-stretch). That this treatment is approximate 
will readily be seen by noting that anharmonic ity terms (cubic, quartic, 
etc.) are omitted from the potential energy. Calculation of precise 
quadratic force constants requires the use of "mechanical" frequencies 
of vibration, the frequencies which would be observed if there were no 
anharmonic terms in the potential energy. This calculation requires a 
substantial quantity of data regarding overtone and combination levels 
of the molecule and was accomplished for HOH vapor originally by Darling 
and Dennison, j^lj 

The four valence force constants were computed here using five fre- 
quencies ( Y ^ and Y 3 of HOH, Y ^ Y 2> and 'Z 3 of DOD) and alternately 
with seven frequencies, (the original five plus Y\ and V 3 for HOD). 

A best fit for each calculation was obtained using a least squares tech- 
nique. The computation was performed on a Control Data Corporation 1604 



TABLE III 



Results of Calculations 
All Frequencies in cm* 1 



Species 


Mode 


Observed 


Calculated 








A 


B 


HOH 


y i 


3614 


3608 


3617 




^3 


3705 


3721 


3705 


DOD 


Vi 


2642 


2645 


2646 




y 2 


1180 


1180 


1180 




^3 


2751 


2725 


27 51 


HOD 


ri 


2692 


2704 


_ __ 




y 3 


3660 


3655 


— 




Force 


Constants* 






Constant 


A 


B 


C 






7.24 


4.67 


8.43 






0.92 


3. 39 


0.77 






0.07 


o 

m 

• 

CM 

1 


o 
1— 1 

o 

1 






-0.68 


-0.84 


0.25 





*Units: Millidynes per Angstrom 



Calculation A: Seven frequencies input (HOH, DOD, and HOD) 

Calculation B: Five frequencies input (HOH and DOD only) 



Calculation C: 



Given by Wilson, Decius and Cross as the mechanical 
frequencies for HOH. Jl 7 j 
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digital computer using FORTRAN 60 input. Details of the calculation are 
presented in the appendix. Results are presented in Table III. 

All of the factors discussed above indicate the conclusion that the 
water molecule when dissolved in a non-polar solvent such as carbon tet- 
rachloride or carbon disulfide interacts very weakly with the solvent. 

In summary these factors are: 

(a) It has been shown that it is consistent with all observations 
to regard the 3804 cm” 1 and 3540 cm” 1 bands as envelopes of rotation 
lines, indicating that the molecules are relatively free to rotate; 

(b) The very low frequency shifts observed upon change from vapor 
to condensed phase indicates that the potential field of the molecule 
is nearly undisturbed by molecular association; 

(c) A calculation based on an unassociated molecule is reasonably 
cons istent ; 

(d) The hygroscopicity of carbon tetrachloride and carbon disul- 
fide is related to the low saturation value of water by a mechanism of 
interstitial trapping which is consistent with minimum molecular inter- 
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APPENDIX 



The quadratic force constants for the potential energy expression 
for the water molecule were calculated by fitting the calculated vibra- 
tional frequencies to the observed frequencies by the method of least 
squares. Calculated vibrational frequencies were obtained through 
solution of the secular equation in the form given by Wilson, Decius 
and Cross, 



I GF - E I =0 

| +* ^ I 

where G, E, and F are matrices, jl 7^ G is the Wilson G matrix given by 



G = 



M l + M 3 



■^2 cos 0 
M ^ cos 0 + m 2 

~/*2 s in 0 -M ^ s in 0 



-Aj sin (? 



-A3 sin 0 



p - [^1 +A 2 +2 ^3 < ' 1-cos ^ 



for non-linear triatomic molecules such as HOH, DOD and HOD. F is 
the matrix of quadratic force constants. For water these constants 
are f r bond stretching; 4< > bond angle bending; f stretching inter- 

action; and f ro ^ , stretch-bend interaction. 



T-i r 



l r 2 



F = 



"r or 



2 L 1 



T 2 *< 



'V* 






JE is the third order identity matrix. The machine solution of the 
equations is quite simple: enter with values for £ and F 9 perform 

the indicated product and then diagonalize this product. The resulting 
eigenvalues are the desired fundamental vibration frequencies. 



The least square fit of data was derived as follows: 
Solve the system of equations: 







0 



where g (Zj) are calculated frequencies, f^ are observed frequencies 
and Zj are the force constants. Differentiating: 

_____ = Q 

^ Z J 




Expanding g^ in terms of Zj: 





k 




^ ®i 

z k 



Az k 




neglecting higher terms so that 




<£ Sj 

^ z j 



= 0 



Let 

£ z k - D ik’ AZfc = A k and i k - g° 

Then 

V [ e i - £ °ik A k] D ij - ° 

V - Xt D ij ■ 0 

1 Lk 

E D k £ i - £ D ji D i* A * ■ 0 

i ik 

Thus 

(d' fc ) . - (d'da). = 0 

where now 
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are column matrices 



(D £ ) . and (D DA) 



(D D) is a square matrix, hopefully with an inverse. Then 



(d'd) -1 D6 - A = 0 



A = (D'D)" 1 DC 



so that 



For the second iteration new force constants (Zj) are 



Z . = Z . , , + A 

jnew jo Id 



where all are column matrices. 
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